Non-integral form of the Gross-Pitaevskii equation for polarized molecules 
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The Gross-Pitaevskii equation for polarized molecules is an integro-differential equation, consequently 
it is complicated for solving. We find a possibility to represent it as a non-integral nonlinear Schrodinger 
equation, but this equation should be coupled with two linear equations describing electric field. These 
two equations are the Maxwell equations. We recapture the dispersion of collective excitations in the three 
dimensional electrically polarized BEC with no evolution of the electric dipole moment directions. We 
trace the contribution of the electric dipole moment. We explicitly consider the contribution of the electric 
dipole moment in the interaction constant for the short-range interaction. We show that the spectrum 
of dipolar BEC reveals no instability at repulsive short-range interaction. Nonlinear excitations are also 
considered. We present dependence of the bright soliton characteristics on the electric dipole moment. 



INTRODUCTION 

Basic properties of collective excitations in a dipolar 
Bose-Einstein condensate (BEC) were studied about ten 
years ago, generalization of the Gross-Pitaevskii equation 
for polarized particles were suggested and used [Ull-llill. 
Since when, evolution of collective excitations in the dipo- 
lar BEC with dipoles being strictly parallel to an external 
field has attracted a lot of attention. Rotons and maxons 
are most fascinating features of the polar BEC spectrum re- 
vealed bv a dilute two-dimensional Bose-condensed dipo- 
lar gas 9], f^, (id\. Instability of the three-dimensional 
dipolar BEC leads to shift of the attention on the quasi two- 
and one-dimensional BEC [11]. Dispersion of excitations 
and other properties of ultracold Fermi and Bose molecules 
has been studied lfl2ll-|[T9h. 

In our previous papers we have considered the spectrum 
of polarized BEC including evolution of electric dipole mo- 
ment direction [20]- [24] using the method of quantum hy- 
drodynamics (QHD) developed in Ref.s |[25ll -l30ll. In this 
paper we consider simpler case. It is usually considered at 
studying of polarized ultracold gases, when all dipoles are 
parallel to an external electric field [Il2ll . lfl6ll . 

Derivation of the QHD equations for electrically polar- 
ized BEC including evolution of the electric dipole direc- 
tion was described briefly in appendixes in Ref.s ll20ll and 
I 2in . Full QHD scheme for polarized particles without the 
short-range interaction was described in Ref. [25]. In this 
paper we describe the derivation of the continuity and Euler 
equations for the fully polarized electric dipole moments, 
which are used here to consider the dispersion of collective 
excitations. Our analysis allows to represent the integro- 
differential Gross-Pitaevskii equation for polarized parti- 
cles in a non-integral form. 

One of the most important points we consider in the 
paper is the explicit form of the potential energy of elec- 
tric dipole interaction. Electric dipoles create the electric 
field and the electric field should satisfy to the fundamen- 
tal Maxwell equations. Thus, we present explicit form 
of the potential energy and trace contribution all terms of 



the potential energy in the dispersion of excitations. We 
present the dispersion of electrically dipolar BEC in com- 
parison with nonpolar BEC reveals the Bogoliubov spec- 
trum, characterized by a linear (phonon) dispersion at low 
momenta, and a quadratic dispersion at large momenta 
(which gives the single-particle dispersion at neglecting the 
third order interaction by the interaction radius term [1291 ). 
We find that the dispersion of electrically polarized three- 
dimensional Bose-Einstein condensate shows no instabil- 
ity. 

This paper is organized as follows. In Sec. II we describe 
a representation of the integral Gross-Pitaevskii equation 
in a non-integral form and present corresponding hydrody- 
namic equations. In Sec. Ill we present the derivation of 
the quantum hydrodynamic equations from first principles 
in the approximation corresponding to the Gross-Pitaevskii 
equation. In Sec. IV we consider the excitation spectrum 
of three-dimensional plane waves. In Sec. V we describe 
the bright soliton propagation in electrically polarized BEC 
in two cases: soliton propagation parallel and perpendicu- 
lar to an external electric field. In Sec. VI brief conclusions 
are presented. 



REPRESENTATION OF THE GP EQUATION AS A 
NON-INTEGRAL NON-LINEAR SCHRODINGER 
EQUATION 

In this paper we are going to discuss different form of 
presentations of the GP equation for electrically polarized 
molecules. It, for instance, can be written as 

ihdMr,t) = l-^A + g\ $(r,t) p 



-d^d^ J dr'G^^{r,r')n{r' ,t)^^r,t), (1) 

where m is the mass of molecule, h is the reduced Planck 
constant, d is the electric dipole moment of molecule, and 
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n{r,t) = ^*{r,t)^{r,t) is the concentration of parti- 
cles. The second term in the right-hand side of equation ([T} 
presents the short-range interaction, and g is the interaction 
constant defined via integral of a short-range interaction 
potential. The last term in this equation describes dipole- 
dipole interaction and contains energy of dipole-dipole in- 
teraction Udd written via the Green function G"^. One can 
see that the GP equation ([Tl is a nonlinear integral equa- 
tion, and the integral term contains a nonlinearity. This 
equation can be presented as a set of two hydrodynamic 
equations M2i1 . the continuity equation and the Euler equa- 
tion. Nevertheless, corresponding Euler equation also con- 
tains the nonlinear integral term (see for example Ref. [l^ 
formula (4.9)). Thus, equation ([Tl is very complicated. Be- 
low, in this section, we will show that this equation can 
be represented in nonintegral form. A nonlinear term ap- 
pears instead of the integral term, and we come to the set 
of equations instead of one equation ([T]). Fortunately, these 
two additional equations are linear equations. Solving of a 
nonintegral equation coupled with two linear equations is 
easier than solving of an integral nonlinear equation. 

Potential of the electric field caused by an electric dipole 
appears as ip = — (dV)(l/r) Isill . consequently electric 
field has the form of E = — V(/?= (dV)V(l/r) or it can 
be rewritten in the tensor form = (i'^V'^V"(l/r). 
Energy of interaction of two dipoles comes as Udd = 
— d2E2i, where 62 is the electric dipole moment of the 
second dipole, and E21 is the electric field caused by 
the first dipole acting on the second dipole. The ex- 
plicit form of dipole-dipole interaction energy is Udd = 
-d^dfV^V"(l/r). It can be rewritten via the notion of 
Green function Udd = -d'^d^G"'^ , where G"'^ =d°'d'^^ 
is the Green function of dipole-dipole interaction. Using 
well-known identity 

_ ^ ^ Srrjr^ ^ ^5"^5(r), (2) 

r r-^ 6 

we can represent energy of dipole-dipole interaction in 
the explicit form. This form contains the term propor- 
tional to the Dirac delta function. At mechanical descrip- 
tion of interaction of two dipoles we can put this term 
away. However, it is very important to keep this term in 
the theory of material fields, such as quantum hydrody- 
namics and the Gross-Pitaevskii equation appearing as a 
part of the quantum hydrodynamic description [|29ll . This 
notice becomes especially important because we integrate 
potential of dipole-dipole interaction in equation ([T]l over 
whole space including point r = r'. We get that elec- 
tric field does not satisfy to the Maxwell equations in the 
absence of a term proportional to the Dirac delta func- 
tion. Choosing a frame of reference with z axes parallel 
to d and neglecting the delta function term we come to 
potential energy Udd.z = —d^{l — 3cos^ 'i?)/r'^. It was 
used for generalization of the GP equation on the system 
of Bose particles having electric dipole moment suggested 
in Ref.s [l]], H, H, it is also discussed in review 01211 . 



This equation does not contain whole electric field that 
leads to mistreating of dipole-dipole interaction. The delta 
function term in dipole-dipole interaction may be included 
in the short range interaction, which is also proportional 
to the Dirac delta function in the pseudopotential descrip- 
tion ll32ll . Explicit consideration of including of the delta- 
functional term in the short-range interaction is presented 
in Appendix. This consideration shows that extracting of 
dipole contribution in the short range interaction confirm 
our results based on consideration of full energy of dipole- 
dipole interaction. Moreover, we are going to show that 
the consideration of full energy of dipole-dipole interac- 
tion gives some benefits. One of them has been mentioned 
already. We have deal with the electric field satisfying to 
the Maxwell equations. Explicitly introducing this elec- 
tric field we can rewrite integro-differential equations (the 
Gross-Pitaevskii equation or corresponding equations of 
quantum hydrodynamics) in non-integral form. We also in- 
clude the Maxwell equations in the considering set of equa- 
tions. Despite the fact that numerical methods of solving 
of the integral Gross-Pitaevskii equation for dipolar BEG 
have been developed, we can admit that solving of the non- 
integral Gross-Pitaevskii equation and Maxwell equations 
is much easier and allows to get analytical results. As ex- 
ample, we present an exact analytical solution for the bright 
soliton in electrically polarized BEG. Moreover, in this pa- 
per the interaction constant of the short-range interaction g 
does not depend on molecule dipole value. 

We can include the energy of dipoles in an external elec- 
tric field —dE^xt in the GP equation ([T]l. Electric dipoles 
create the electric field, including their internal field, we 
can rewrite the last term in equation ([T) in form of — dE 
introducing the electric field caused by electric dipoles. In 
result we find 
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A + 5 1 $(r,0 1^ -dE U.(r,t), 



(3) 

where E = ^ext + E,;„t, Ei„f is the electric field caused by 
dipoles. However introducing additional physical variable 
we have to present equation for this variable. Looking on 
the explicit form of internal electric field 



J dr'G"^{r,r')n{r',t) 



we can obtain that this electric field satisfies to the Maxwell 
equations 

VE(r,t) = -47rVP(r,t) 



-47r(dV)n(r,t) = -47r(dV)($*(r, t)$(r, t)), (4) 



and 



VxE(r,t) = 0. (5) 

Equation Q is also called the Poisson equation with the 
charge equal to zero and explicit writing of the electric dis- 
placement field via the electric field E and the polarization 
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P. Equation (|5]) is the quasi static form of the Faraday's 
law of induction. For getting Maxwell equations we have 
to use the full Green function of dipole-dipole interaction. 
Equation ^ is the non-integral form of GP equation for 
electrically polarized particles. 

Non-linear Schrodinger equation (|3} can be represented 
as a set of two hydrodynamic equation, they are the conti- 
nuity equation 

dtn + V(nv) = 0, (6) 
and the Euler equation 

4m 4m V n 



-gnd^n + ndd^ilE) 



(7) 



where 1 is the unit vector in the direction of the polarization 
formed by the external field, n is the particle concentration, 
n = and v is the velocity field. 

In the last term of the equation (Q, ndl is the po- 
larization of medium for fully polarized dipoles with no 
evolution of dipole direction. This term appears instead 
of P^{r,t)d''E^{r,t), where P is the polarization of 
medium for general case including the dipole direction evo- 
lution around the equilibrium position formed by the exter- 
nal electric field. Below we present the derivation of equa- 
tion (|7]l from a microscopic theory (from the many particle 
Schrodinger equation). 



DERIVATION OF QHD EQUATIONS FROM 
SCHRODINGER EQUATION 

We have discussed the Gross-Pitaevskii equation for po- 
lar molecules and corresponding equations of hydrody- 
namics as a master equations. In this section the deriva- 
tion of the quantum hydrodynamic equations from the first 
principles is described. 

Starting from the many-particle Schrodinger equation 

tndtip{R,t) = Htl;{R,t), 
(ri, r^v) is the set of coordinate of N parti- 



where R 
cles, with the Hamiltonian 



H 
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(8) 



we construct a system of QHD equations for particles hav- 
ing an electric dipole moment . The following designa- 
tions are used in the Hamiltonian (HJ: = —ihd", E"'^^^ 
is the electric field, m^ is mass of particles, h is the reduced 
Planck constant, and Gff = d°'d^ 1 /vij is the Green func- 
tion of dipole-dipole interaction, Uij = U{\ri 



potential of the short-range interaction. Writing Hamilto- 
nian dH we suppose that particles have electric dipole mo- 
ment, but directions of electric dipole moments are, in gen- 
eral case, different, but they can be polarized in the direc- 
tion of external field. The dipole-dipole interaction leads to 
the evolution of both positions of particles and directions of 
the electric dipole moments. 

The first step in the construction of the QHD method is 
to determine the concentration of particles in the vicinity 
of r in the physical space. If we define the concentration 
of particles as quantum average of the concentration oper- 
ator in the coordinate representation h = ^ . 5{r — r^) we 
obtain 



n 



:r,t) = j dRY,Kr-r,)r{R,t)'^{R.t), (9) 



where dR = H^i ^Vp. 

Differentiation of n(r, t) with respect to time and apply- 
ing of the Schrodinger equation with Hamiltonian (|8) leads 
to the continuity equation + Vj = 0, where the current 
density takes the form of 

j(r,t) = j dRY,5{v-v;)^{ij*{R,t){vMR,t))+h.c. 

(10) 

In hydrodynamics one usually has deal with velocity field 
v(r, t) instead of the current density j(r, t), thus we have 
to present a connection of j and v. 

The velocity of i-th particle Vi(i?, t) is determined by 
equation 
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v,(i?,t) = —V,S{R,t 



(11) 



where t) presents the phase of the wave function 

'iS{R,ty 



4^{R, t) = a{R, t) exp 



Tj I ) is the 



The quantity Vj(i?, t) describes the current of probability 
connected with the motion of i-th particle, in general case 
yi{R, t) depends on coordinates of all particles of the sys- 
tem R, where R is the totality of 3N coordinates of N 
particles of the system R = (ri, rjv). Velocity field 
v(r, t) is the velocity of the local centre of mass and deter- 
mined by equation j(r, = n(r, t)v(r, t). This means that 
Ui(r, R, t) = yi{R, t) — v(r, t) is a quantum equivalent of 
the thermal velocity. 

A momentum balance equation can be derived by differ- 
entiating current density ([TOt with respect to time: 

dtr{r,t) + -d^U"^{r,t) = -F"(r,t), (12) 
m m 

where F" (r, t) is a force field and 11"'^ (r, t) is the momen- 
tum current density tensor. 

Performing explicit separation of particles' thermal 
movement with velocities Ui{r,R,t) and the collective 
movement of particles with velocity v(r, t) in equations of 
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continuity and of the momentum balance (fT2] l we can find 
that the tensor 11"^ (r, t) takes form of 



n"^(r,t) 



mnv 



In this formula p"^(r, t) is the tensor of kinetic pressure. 
This tensor tends to zero by letting — > 0. The tensor 
T"*^ is proportional to H"^ and has a purely quantum origin. 
For the system of noninteracting particles, this tensor is 

T"'^ = -^d^d^n + ^-(d"n)(d^n). (13) 
4m 4m n 

This term is called the quantum Bohm potential. 

As the particles of the system under consideration in- 
teract via long-range forces the approximation of the self- 
consistent field is sufficient to analyze collective processes. 
With the use of this approximation two-particle functions 
appiaring in the momentum balance equation can be split 
into a product of single -particle functions. Taken in the ap- 
proximation of self-consistent field, the set of QHD equa- 
tion, continuity equation and momentum balance equation 
has the form 



and 



mn 



dtu + V(nv) = 0, 



{dt + vV)7;" + dpT''^ = F^^ 



(14) 



SRI 



+ P^d"Eit + P^d'' I dr'G^^r,r')P^r',t). (15) 

Let's discuss the physical meaning of terms on the right- 
hand side of (ITSl ). Force density of the short-range interac- 
tion i^5fl/(r, t) was considered in Ref. ll29ll . It was found 
that in the first order by the interaction radius the force 
density has form corresponding to the Gross-Pitaevskii ap- 
proximation Fgjjj{r,t) = —gnd^n. The second term is 
the effect of non-uniform external electric field on the po- 
larization density. It should be noted that the form of this 
term is distinct from the expression that describes the force 
affecting a single dipole. The third term of the equation 
( fTsl l describes a force field that represents interactions be- 
tween particles, namely the dipole-dipole interaction. 

Note, that for a three-dimensional system of parti- 
cles the momentum balance equation ( [TSl l may be writ- 
ten down in terms of electrical intensity of the field 
that is created by dipole moments of the particle sys- 
tem: P^{r,t)d°'E^{r,t), where E™t(r,t) satisfy to the 
Maxwell equations Q and (|5]l. 

The method we develop in this work is valid both for 
bosons and fermions. The type of statistics that particles 
are subject to affects the calculation of many-particle func- 
tions (correlations) that evolve in the momentum balance 
equation dTSl l and are neglected in the self-consistent field 
approximation. A method for the calculation of correla- 
tions in the QHD equations has been developed in papers 



The polarization evolves in the momentum balance 
equation ( fTSl) is 



P"(r,t)= / dRy6{r-r,)r{R,t)dtiP{R,t). (16) 



This function presents general definition of the polariza- 
tion appearing in the quantum hydrodynamics. It can be 
replaced with ndl at assumption that all dipoles are paral- 
lel to each other. If we need to consider the dipole direc- 
tion evolution we can derive equations for the polarization 
P, and, therefore, expand the set of the quantum hydrody- 
namic equations (see for example [201 . 123.1 . ll25ll ). 



EXCITATIONS 

calculation of the excitation dispersion 

Solving the set of QHD equations 0, Q, (l6|, and © at 
assuming that all dipoles are parallel to the external elec- 
tric field, and the inter-particle interaction does not change 
the direction of dipoles, we consider wave perturbations 
around the equilibrium state described with the following 
parameters n = Uq, v = 0, E = E^xt- In the lin- 
ear approximation on small simple harmonic perturbations 
6f = f{u},k)exp{—iujt + ikr) of the particle concen- 
tration 5n, the velocity field 5y, and the electric field 6E, 
where w is the frequency of the wave oscillation and k 
is the wave vector. Amplitudes of the wave perturbation 
satisfy to an algebraic form of the quantum hydrodynamic 
equations: the continuity equation ^ 

— iu}6n{uj, k) + noik6\{uj, k) = 0, (17) 

the momentum balance equation d?]) (the Euler equation) 

—iujmnQ6v°'(uj,k) H zk"'k^6n(iLi,k) 

4m 

= -gnQik'^6n{oj, k) + mQdk'^\5E{oj , k), (18) 
and the Maxwell equations and ^ 

kxEx{u),k)+kyEy{uj,k)+kzEz{oj,k) = —A'Kdkz5n{(jj,k), 



(19) 



and 



k X (5E(a;,k) = 0. (20) 

The continuity equation ( fTTl) gives a connection of the 
particle concentration amplitude 5n{u),k) with the ampli- 
tude of the velocity field (^v(w, k): 

k(5v((jj, k) 



5n{ijj, k) = tiq- 



Using this relation we find connection of the velocity field 
amplitude with the amplitude of the electric field from the 
Euler equation ([Tsl l, which is 



(k(5v(w,k)) 



mui 



-u}dP6E,{ u},k) 
gnok'^ 



2 _ ft^fc" 



(21) 
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We will put this formula in equation ( [T9l l, but now we ex- 
press 5Ex and SEy via 6Ez using equation (l20l) . We have 



6E, = ^6E,, 



and 



5E„ 



'-5E^. 



(22) 



(23) 



These formulas we also put in the formula ( [T9l l. The for- 
mulas (|22] | and ( |23l ) are uncorrect if = 0. If A;^ = we 
have E^ = 0. Using it we find the Bogoliubov spectrum 
because electric dipoles give no contribution in the Euler 
equation, and consequently, in the particle evolution. 
In result the equation ([T9l l appears in the following form 



k^5E, 



And^klk^nodE, 



4m 



(24) 



kl + kl + kl 



where k'^ 

Assuming that the amplitude of the electric field is not 
zero we cancel SE^ and get the dispersion equation 



„ h^k^ qrink'^ AirnnCpk'i 
= — ^ + + -. 



(25) 



This formula gives the dispersion of collective excitations 
in the electrically polarized three-dimensional BEC at as- 
sumption that all dipoles are parallel to the external electric 
field and the dipole-dipole interaction causes no evolution 
of the electric dipole direction (the polarization waves is 
the example of the coherent evolution of the dipole direc- 
tion il, iH, iQ). 



properties of the excitation dispersion 

Our calculation gives the following spectrum of the col 
lective excitations 

h^k^ griok'^ A-KUi^dPk'^ cos^ 6 



UJ 



+ 



+ 



Am? m 

where cos 9 = k^/k. 

This result differs from formula 



m 



(26) 



obtained in the previous papers (see, for example, formula 
(5.1) in Ref. yjj or formula (11) in Ref. lil6J), is 
the dipolar coupling constant C^ii = d^/eo in the SI units, 
where is the vacuum permittivity, or C^^i = And^ in 
the CGS units. In formula (IZTI l the constant of the short- 
range interaction depends on the electric dipole moment of 
molecule g = g(d) [16] (see formula (3) and the text after 
it). To get correct properties of the spectrum we need to ac- 
count this dependence explicitly (we present this consider- 
ation in the Appendix), but these properties have been con- 
sidered with no account of dependence g{d). We present 



the brief discussion of formula (|27] |. For unpolarized BEC 
g > corresponds to the repulsive short-range interac- 
tion. The attractive short-range interaction g < leads 
to the unstable spectrum of unpolarized three-dimensional 
BEC. Let us discuss changes of the dispersion in the re- 
pulsive BEC g > caused by dipoles. The contribution 
of dipoles in the dispersion dependence may have differ- 
ent sign depending on 9. As follows from formula (|27^ 
square of frequency become negative at large enough Cdd 
and some angles, that reveals an instability of the three- 
dimensional polarized BEC. Under proper conditions the 
spectrum (l27l l reveals a minimum in the dispersion at in- 
termediate momentum |6]. The formula (|27^ is obtained 
using the potential of dipole-dipole interaction in the ab- 
sence of 5-function term. 



Let us discuss the spectrum (|26] | obtained in this paper. 
At small wave vectors k the frequency can be written as 
LO = {nQ/m){g + Aird'^ cos^ 9)k. We can see that, at 
the repulsive short-range interaction g > 0, the coefficient 
g + 4:7rd^ cos^ 9 is positive for all angles 9. If And^ /g > 1, 
at small angles 9 <ti tt/2 the quantity g + And^ cos^ 9 is 
positive even for the repulsive short-range interaction g < 
0. It shows us that the electric polarization of BEC does 
not lead to destabilization of the spectrum of BEC. 



BRIGHT SOLITON IN THE ELECTRICALLY 
POLARIZED BEC 



The fundamental nonlinear excitation in attractive (g < 
0) unpolarized BEC is the bright soliton, which is an area 
of the increased concentration. We consider a perturba- 
tion of the hydrodynamic variables n, v, and E as functions 
of X — Vot and assume that external electric field applied 
along z axes. Therefore we consider propagation of the 
bright soliton perpendicular to the external field. We find 
that electric field gives no contribution in evolution of par- 
ticle concentration n and velocity field v. The last term in 
equation ^ is equal to zero, since dxE^ = that follows 
from equation dU. Consequently, we have no change in 
properties of the bright soliton. 

Next, we consider the bright soliton propagation parallel 
to the external field, so we suppose that n, v, and E are 
functions of ^ = z — VqI. In this case the last term in the 
equation (|7]) appears as ndd-^E^ = —ATrd^ndzn that we 
find using equation (O. Thus we can combine two terms 
in the right-hand side of the equation (|7]l, which present z 
projection of the force density F, and write it in the form 
of ^2 = —{g + 4Trd^)ndzn. We have that the dipole- 
dipole interaction of fully polarized electric dipoles leads 
to change the interaction constant. All difference, we find 
in compare with unpolarized case, is replacement of g with 
g = g + Aird"^. Using the methods described in Ref.s (3^, 
1 34j] we come to the particle concentration profile in well- 
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known form of 



n{0 = no 



1 



cosh a 



(28) 



where Uq = 2m | | /(| ^ |), a = a/2 | E \m£^/h, 
where E = Cq — Vq/2 and Vq is the velocity of the soliton 
propagation, Cq is a constant. The bright soliton exists at 
^ < and 5 < 0. 

The bright soliton exist at the condition g + AirdP < 0, 
which shows that the polarization leads to decreasing of the 
area of the bright soliton existing. In the electrically polar- 
ized BEC module of the short-range interaction constant 
must be more than AttcP. 



CONCLUSION 

Explicitly considering the correct potential energy of 
the electric dipole-dipole interaction and assuming that all 
dipoles are parallel to the external field we presented a non- 
integral form of the GP equation for electrically polarized 
molecules and corresponding hydrodynamic equations (the 
continuity equation and the Euler equation). Solving these 
equations we obtained spectrum of collective excitations 
in the electrically polarized BEC. We found that the ob- 
tained spectrum differs from the one presented in litera- 
ture. 3D uniform dipolar BEC does not reveal any insta- 
bility for all directions of wave propagation. Introduction 
of the internal electric field satisfying the Maxwell equa- 
tions allows to represent the GP equation in a nonintegral 
form and gives possibility to study the bright soliton an- 
alytically. The bright soliton in the electrically polarized 
BEC was considered using described approximation. We 
shown that area of parameters of the bright soliton exist- 
ing is decreased in compare with the unpolarized BEC, but 
its form does not change, for the soliton propagation along 
the external field. And the polarization gives no influence 
on the bright soliton properties at the soliton propagation 
perpendicular to the external field. 



APPENDIX 

Writing "Schrodinger equation" for the condensate wave 
function, which appears as a part of ip operator ll32ll . we 
have 



+ K.t(r) + y $*(r',t)C/(r-r')$(r',t)dr'j$(r,t) 



Particles interact by means a short-range interaction and 
the long range dipole-dipole interaction, then 

U = g6{r-r') + —d'5{r-r') + UM„ (30) 



where g = iirh^a/m, a is the scattering length describing 
short-range interaction existing at c? = 0, and 

dd-3(rd)(rd) 



U, 



ddO 



(31) 



Delta functional terms can be combined U = g'd(r — r') + 
Uddo, where g' = g'{d) = g + 47rd^/3. If amplitude 
of dipole-dipole interaction proportional to dP is large in 
comparison with the short-range interaction g <^ d^, when 
we can write g' = And^ /3. Using explicit form of g'{d) 
and Cdd we rewrite formula (ETt and find 

w = [g + —d + — — (3cos 61-1 +——, 

(32) 

combining together terms proportional to d^ we get for- 
mula ( l26l) obtained in this paper 

The author wish to thank Professor L.S. Kuz'menkov for 
discussion of the results obtained. 
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